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Abstract. In the recent papers [i, 0] the authors study the existence of weak 
solutions and the hydrodynamic Umit of kinetic flocking equations with strong 
local alignment. The introduction of a strong local alignment term to model 
flocking behavior was formally motivated in these papers as a limiting case 
of an alignment term proposed by Motsch and Tadmor [ ]. In this paper, we 
O rigorously justify this limit, and show that the equation considered in [1, -5] 

is indeed a limit of the Motsch- Tadmor model when the radius of interaction 
goes to zero. The analysis involves velocity averaging lemmas and several 
estimates. 

o 



< 



1. Introduction 

PLh In [G] Motsch and Tadmor identify an undesirable feature of the widely studied 

Cucker-Smale flocking model (cf. [1, 2, 3]): In the Cucker-Smale model, the align- 
ment of each individual is scaled with the total mass such that the effect of alignment 
is almost negligible in sparsely populated regions. To avoid this, they propose a 
new model in which the alignment term is normalized with a local average density 
y instead of the total mass. Motivated by this work, the authors of the present paper 

proposed in [4] to combine the Cucker-Smale and Motsch- Tadmor models, letting 
the usual Cucker-Smale alignment term dominate the large scale dynamics and the 
^ Motsch- Tadmor term the small scale dynamics. This remedies the aforementioned 

deficiency while maintaining the large scale dynamics of the Cucker-Smale model. 
At the mesoscopic level, the proposed model takes the following form 

ft + div,(w/) - div,(/V,*) + div,(/F[/]) + div„(/L'-[/]) = 0. (1.1) 

Here, the unknown is the distribution function / := f{t,x,v). The first alignment 
term F[-] is the standard Cucker-Smale alignment term given by 

> F[f{x,v)]^ f <P{x-y)f{y,w){w~v) dwdy, (1.2) 

where $(a;) is the influence function (e.g ^{x) = l/{l + \x\'^)). The second alignment 
^ term L''[-] is the Motsch- Tadmor alignment term given by (see [ii]): 

= /g^d K'^jx - y)f{y, w){w - v) dwdy 
jji2aK'-{x-y)f{y,w) dwdy 
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where the index r denotes the radius of influence of (see (1.4) below for the 
definition of K^). Finally, the function '^{x) is a given confinement potential intro- 
duced to avoid mass vanishing to infinity (it satisfies lim|2,|_j.oo ^'(x) = oo). Note 
that this term is not necessary if we assume, for instance, that the initial distri- 
bution /(O, a;,u) has compact support in x and v (since (1.1) will propagate this 
property) . 

The only difference between (1.2) and (1.3) is the renormalization by the local 
average density J^2d K^{x — y)f{y, w) dwdy. We can also write U as follows: 

where 

^r^^^ ^ /R2d K^'jx - y)wf{y, w) dwdy 
/R2d K'-ix - y)f{y, w) dwdy ' 

In this form, it is obvious that the strength of the alignment force is now indepen- 
dent of the total mass, which was the original intend of [ ]. Another effect of this 
renormalization is to break the symmetry of the alignment. As a consequence, (1.1) 
does not conserve momentum nor energy, and the derivation of an energy bound 
will be one of the main difficulty in the analysis of (1.1). 

The purpose of this paper is to study the limit r ^ in Equation (1.1) when 
the function converges to the Dirac distribution 6o. In other words, we study 
the limit of (1.1) when the Motsch-Tadmor term div^,(/L''[/]) becomes a local (in 
space) alignment term. For the sake of simplicity, we assume that if has the form 

K'-ix) ^r-'^K (^^y (1.4) 

where if is a given function satisfying 

K e Cc(M''), K{0) > 0, / K^'ix) dx = 1. (1.5) 

Js.'' 

When r — > 0, we then formally expect to have 

~rr -.r^O , . lM^wf{x,w) dw 

u [x] — > u{x) — 



jRdf{x,w) dw 



and so 



J^df[x,w) dw 

Passing to the limit in (1.1), we thus obtain the equation 

ft + div,(i;/) - div„(/V,*) + div,(/i^[/]) + div„(/(u - v)) = 0, (1.7) 

which is studied in [ , ]. The new local alignment term can also be seen as a local 
friction term centered at u. 

The purpose of this paper is to rigorously justify this limit r — > 0. More precisely, 
we will prove the following theorem: 

Theorem 1.1. Let < fo G ^^(K^d) n L°°(M2d) ^j„g„ ^nd T be a finite final 
time. For each r > 0, let f^(t,x,v) be a weak solution of (1.1) in the sense that 

I -r<Pt-vf'V^<P + f'V^-^V^<P-rF[f']V,<P dvdxdt (1.8) 
- / /'■^'■[/nv„0 dvdxdt = / /o?!)(0, •) dvdx, V</) e C^{[Q,T) x W"^), 
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where U" is given by (1.3) and is given by (1-4). Then, as r — > 0, 

inL°°(0,r;L°°(M2'i)nLi(E2'i)), 

rL^[n f{u - v) m L%{0, T) x M^^)^ g < 

d+1 

with u(t,x) such that j — pu (see (2.1) for a precise definition). Furthermore, the 
limit f(t,x,v) is a weak solution of (1.7) in the sense that 

I ~f(t>t - vfV^(t> + fV^^V^<t> - fF[f]V^<t) dvdxdt (1.9) 
f{u - v)V„0 dvdxdt = / /V(0, •) dvdx, \f(j) e C^{[0,T) x M^d), 



2. Preliminary material 

In this section we have gathered some results that will be needed to prove The- 
orem 1.1. We begin by introducing some convenient notations. We denote the first 
and second moments of /, and their weighted counterparts, as follows: 

Q{t,x)^l f{y,w) dwdy, Q''{t,x)^ K'' (x - y)f{y,w) dwdy, 

j{t,x)= f{y,w)w dydw, f{t,x)= K'' {x - y)f{y,w)w dwdy. 
We also define the corresponding velocities 

U(t,X) — — r, U [t,X) — r. 

g[t,x) Q^(t,x) 

Note that the definition of u (and u) is ambiguous if p (resp. g) vanishes. We thus 
define u pointwise by 

. ^ f 4^ lip{x,t)^Q , , 

u{x,t) = \ p{x,t) , 7^ (2.1) 

[ if p(a;,i) = 

Since we have 

1/2 



J<(^J\v\^fix,V,t)dv^ 



the bound on the energy of / will imply that j — whenever p — and so (2.1) 
implies in particular j = pu. 

With the above notation, we have L'^[f] —u^~v, and (1.1) can be written as 

ft + div,(/t;) - div„(/V,*) + div,(/^^[/]) + div,(/(2'- - v)) = 0. (2.2) 

The following proposition states that (2.2) is well-posed in the sense of weak 
solutions (see [ ] for the proof). 

Proposition 2.1. Assume that < /o G [L°° n L'^]{R'^'^) and T < +oo are given. 
Then, for any r > 0, (2.2) admits a weak solution < / G C{0,T; L^(R'^'^)). 
Moreover, f satisfies 

CT 

ll/llL°=(0,T;LP(R2<i)) < e P' |!/o|!LP(R2<i), (2.3) 

[ {\v\^ + '^ix)) fit,x,v)dvdx <Ce^^£iO), (2.4) 

where the constant C might depend on r. 

To conclude this section, we recall the following classical lemma, which will be 
used to derive L'p integrability of g and j (see [ ] for the proof): 
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(2.5) 



Lemma 2.2. Assume that f satisfies 

ll/llL-([o,T]xR2<i) < A^, and / \v\'^ f dvdx < M. 

Then there exists a constant C = C{M) such that 

IIpIIl-(o,T;Lp(R'')) < C, for every pe [1,^), 
IIjIIl~(o,T;Lp(R'')) < for every p e [1, ftf ), 

where p — J f dv and j — J vf dv. 

2.1. The Velocity Averaging Lemma. When passing to the limit in (2.2), the 
main obstacle is to obtain compactness of the product fv/". The instrument we 
will use to obtain this is the celebrated velocity averaging lemma. We will use the 
following version due to Perthame & Souganidis [ ]. 

Proposition 2.3. Let {/"}n be bounded in L^^^{R'^'^'^^) with 1 < p < oo, and 
{G"}„ be bounded m L[„^(M2'i+i). ///« and G" satisfy 



loc\ 

,^^)2(i+l^ Tf fn r<ri 
loc\ 



/r + V ■ V,/" = v:;g", r\t=o = /" e (i 

for some multi-index k and ip e Clf^^R'^'^), then {g^^} is relatively compact in 

The previous proposition cannot be directly applied to obtain the needed com- 
pactness. In fact, we will rely on the following lemma which can be seen as a 
corollary of the previous proposition. The proof can be found in [4]. 

Lemma 2.4. Let {/"}« and {G"}„ be as in Proposition 2.3 and assume that 

/" is bounded m L°"{R^'^+^), 

(|i'P + ^')/" is bounded in L°°(0, T; ^^(M^'i+i)). 

Then, for any (p{v) such that Ifiv)] < c\v\ and q < the sequence 

r^iv)dv\ , (2.6) 



is relatively compact in L^((0,T) x M.'^). 

2.2. An important technical lemma. In view of Lemma 2.2 and 2.4, it is clear 
that in order to get convergence results for f^ and its moments, we will need to 
obtain some estimate on f^ that are uniform with respect to r. The main difficulty 
will be to show that the energy estimate (2.4) holds with constants independent on 
r (which does not obviously follows from the result of [4]). For this we will make 
use of the following technical lemma, which can be found in [ '] (the proof is given 
below for completeness): 

Lemma 2.5. Assume that there exists < i?i < -R2 < oo such that 

K{x)>0 for\x\<Ri, K{x)=0 for\x\>R2. (2.7) 
There exists a constant 

such that 

for all nonnegative functions p G L^{M.'^). 
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The most important part of this lemma is the formula (2.8), which implies that 
if we replace the function K with aK{f3x), for any a > and /? > 0, then the same 
results holds with the same constant. 

We deduce: 

Corollary 2.6. Assume that is given by (1.4) where R satisfies (1.5). Then, 
there exists a constant C independent of r such that 

^ ''^ kiK'-{\x~ z\)p{z)dz'^^ VyeM 
for all nonnegative functions p G iy"'^(M''). 

Proof of Lemma 2.5. We recall that p{x) — Jj^^ K{x — z)p{z) dz and we note that 



/ K{x -y)^dx< (sup K) f 



dx. 



,Ay) Pi^) 

Next, we cover Bji^iy) with balls of radius Ri/2: We have 

N 

BbM C \JBr^/2{x,) 
i=l 

with N ^ {R2/R1Y. We can thus write 

/ K{x-y)P^d:,<{su^K)j2l ^-^,dx. 
Moreover, clearly, 

p{x) = / K{x — z)p{z) dz > K{x — z)p{z) dz. 

Jr" Jsa^/^ix,) 

By combining the two previous inequalities, we see that 

N 



[ K{x-y)^dx<{snpK)Y, f 



K{x — z)p{z) dz 



dx. 

Now, using the fact that when a;, z e B}i^/2{xi) we have \x — z\ < Ri, we deduce 



/ K{x~y)'Jfld:r<^-P^Y.f 1 '-^^^^^ dx 



supi^r ^ supiC / R 



J p{z) dz 
d 



inf Bflj (0) K inf (0) K \Ri ^ 

and the proof is complete. □ 

2.3. A priori estimate. We can now conclude this preliminary section by proving 
that satisfies some a priori estimates uniformly with respect to r. We recall that 
the energy functional is defined 

f W = / f ^ + ^(a^)') /(i, 2;, v) dvdx. (2.9) 

We then prove: 

Proposition 2.7 (Energy bound) . Let < /o G {M?'^) r\ {M?'^) he given, letT 
be a finite final time, and let f be the corresponding weak solution of (1.1). There 
is a constant C > independent of r > such that .such that 

CT 

ll/llL°=(0,T;LP(R2d)) < e^||/o|lLP(R2d), (2.10) 
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and 

sup £{t) + l^ I f\u~v\'^dvdx 
te(o,r) 2 7R2d 

1 /■ V ■ 

+ o / "^(^ - - ""1^ dwdydvdx < C{T)£{0). 

The proof of Proposition 2.7 relies on two auxiliary results (Lemmas 2.8 and 2.9 
below) which we will prove prior to proving the proposition. We begin with the U' 
estimate (2.10): 

Lemma 2.8. Let f be a weak solution of (1.1). There is a constant C , independent 
of r, such that 

sup \\f\\Lv(K?i) < ||/o||LP(R2<i)e^^. (2.12) 
te(o,T) 

Proof. Let B{f) be a continuous function and let b{f) = fB'{f) — B{f). By 
multiplying (1.1) with B'{f) and integrating, we obtain 



d_ 
dt 



[ B{f) dvdx = I v\/^b{f) dvdx 

jR2d J^2d 

+ f {F{f) + U{f)-V,-^)VMf) dvdx (2.13) 
b{f) (div,„ F{f) + div„ U{f)) dvdx. 
Next, using the definition of the alignment terms, we see that 
divi, F[f] ^ -d ^{x - y)f{y, w) dwdy, 

div = ^ /R2d K"'{x-y)f{y,w) dwdy ^ 

/R2d K'-ix - y)f{y, w) dwdy 

Setting these identities in (2.13), we find that 

d 
dt 



B{f) dvdx ^ / b{f)ld + d ^{x - y)f{y,w) dwdy ] dvdx 
6(/)((i + dM||<f>||ioc.(R2d)) dvdx, 



< 



where Af is the total mass. Next, we let B{f) = fP such that b{f) = (p— 1)/^- An 
application of the Gronwall inequality then provides the bound 



sup ||/|lLP(R2d) < II /oil LP (R2d)e p 
te(o,T) 



-CT 



which is what we set out to prove. 

□ 

The main difficulty in proving the energy estimate (2.11) (even for r > 0) is to 
control the non-symmetric Motsch-Tadmor alignment term. This is the goal of the 
following Lemma, which relies on Lemma 2.5: 

Lemma 2.9. There is a constant C, independent of r, such that 

[ fvU'[f] dvdx < C£{t) - i / f\u' - dvdx. (2.14) 

jR2d 2 jR2d 
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Proof. By definition of L'' , we have tliat 

L''[/] = — / K''{x-y)f{y,w){w-v) dwdy 

'^l^ . (2.15) 

K''{x - y)[j{y) - g{y)v) dy ^ L - v:=u - v. 



By adding and subtracting, we obtain 

/ fvL^lf] dvdx — / f(u — v)v dvdx 

= -l I f[u-vf dvdx+\ [ fu^ - fv^ dvdx (2.16) 
< — - / f{u~v)'^ dvdx + - j gu^ dx. 
From the Holder inequality, we have that 

eu-= K'^{x-y)f{y,v)v dvdy <gi [ / K''{x - y)f{y,v)v'^ dvdy 

■/K^d 

Hence, the following inequality holds 

gu^ < / K^{x — y)f{y,v)v'^ dvdy, 

from which we deduce 

gu^ dx < K^{x — y) f{y,v)v'^ dydvdx 

g{-) (2.17) 



< sup f / K'-ix ~y)^ dx] E(t) < C£{t), 
V kJr g[x) ) 



where the last inequality follows from Lemma 2.5. Inserting (2.17) in (2.16) con- 
cludes the proof. 

□ 

We have now gathered all the ingredients we need to prove Proposition 2.7. 

Proof of Proposition 2.7. Only (2.11) remains to be proved. By direct calcula- 
tion, 

f^ + f^-^dvdx= ft<f + ft^-^ dvdx (2.18) 
at 7R2d z Jp2d z 

- / vfV.^iJ - vfV,^ + fvF[f] + fvU[f] dvdx. 
Using the symmetry of we write 

fvF[f] dvdx — I $(x — y)f{x, v)f{y, w){w — v)v dwdydvdx 



$(x — y)f{x, v)f{y, 'w){v — w)w dwdydvdx 
$(x — y)f{x, v)f{y, w)\w — up dwdydvdx. 
Then, we conclude the proof by applying this identity and Lemma 2.9 to (2.18). 



□ 
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3. Convergence and proof of Theorem 1.1 

Equipped with the bounds of the previous section, we are ready to send r — > in 
(1.1) and thereby proving Theorem 1.1. For this purpose, we let {?""}„ be a sequence 
of positive numbers such that r" — > as n ^ cxi and consider the corresponding 
solutions /" of 

/r + div,(«r) - div„(rv,vi;) + div,(rF[/"]) + dw^^iu" - v)) = o, (s.i) 

where we recall the notation 

_ /R2d K''"'{x - y)P{y,w)w dwdy 
' /R2d K''"{x - y)r'{y, w) dydw 

Our starting point is that Lemma 2.8, Proposition 2.7, together with Lemma 2.2, 
asserts the existence of a fimction < / G C(0, T; L^(M2d)) p x°°(0, T: L°°(M2rf))^ 
such that, as n — oo, 

in L°°{0,T;L°°{R^'^)nL\R^'^)), 

d + 2" 



u = 



^ g in L°°((0,r);LP(M'')), for every p £ 
j" ^ j in L°°{{0, T); LP(M'*)), for every p e 



d J' (3.2) 

d + 2" 



d+1 



Moreover, the velocity averaging Lemma 2.4 is applicable. By setting ip{v) — 1 and 
^p{v) = u in Lemma 2.4 we obtain respectively 

g"^g in L«((0,T) X M'^), for every g < 
j" ^ j in L«((0, T) X R'^), for every q<^^^ 



(3.3) 



d+V 

along some subsequence as n — )■ oo. Furthermore, we can prove: 

Lemma 3.1. Given the convergences (3.2) - (3.3), we have 

~r^J, r^g, L'((0,r) xM'^), (3.4) 

where the convergence takes place along the same subsequence as in (3.3). 

Proof. We commence by recalling the following classical results concerning moUi- 
fiers like = : For any e > 0, there is a m such that 

\\K" -k g- gWLiiR") < e, n>m. 

Now, consider a subsequence n*^, where > m, along which — > q. By adding 
and subtracting, we obtain 

11^ - (?llL5(Krf) < lie*" - i^"*£'||^,(Rd) + £'|lL<j(Rrf) 

= ||i^"*(£'" - e')llL9(Rd) + WK"^* g~ g\\Li(K.'i) 

< l|(?"-£'llL.j(R^)+e = 2e, 

for any q < The same argument can be applied to prove compactness of j", 

which concludes the proof. 

□ 

Lemma 3.2. From the convergences (3.2) - (3.3), it follows that 

d+2 



/"y"A/y tn L°°{{0,T);LP{R'^)) for every pG 



^'d+1 
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Proof. For the sake of clarity, let us introduce the notation 

For any smooth function 'ijj{x,v) := (f>{x)ip{v), we write 

dvdx = I u"0(a;) ( / dv ] dx 



(3.5) 

u^g^cf) dx = rnJ^cj) dx. 
Now, using the Holder inequality, we find that 

Wm'^WLHR-^) < II¥'IIl-(r'')II(?"P^, ,JI((?")^w"I1l2(r'') 

<C||e"||t(^,)||(f.")5S"|U.(K.), 

which is bounded by (2.17) and Lemma 2.2 provided 

d+2 d+2 

Hence, there exists a function m e i°°((0, T); L''(M'')) and a subsequence such that 

m^^m in L°°((0,r);LP(M'')), for every p € 1.^^^ 

and it only remains to prove that 

m — ug^p, where u is such that j = gu. 
Let us first verify the existence of such a function u. Consider the set 
Ar = {{t, x) e Sfi(O) X (0, T); g{t, a:) = 0} , 
where -B_r(0) is the ball of radius R centered at 0. By direct calculation, 

/ \in\ dxdt < ( [ dxdt] [ [ g" dxdt] 

<Ct(1 



and hence we have that j = a.e in A^. If we define the function u as 

|0, iig{t,x) = Q, 

we have that j — gu and it remains to prove that m — g^u. To this aim, we first 
observe that we can deduce as in (3.6) that 

WKWmA.) <c\mhiA.)^^ 

and hence it suffices to check that 

m{t, x) — u{t, x)g^p{t, a:), whenever g{t, x) ^ 0. 

For this purpose, we consider the set 

B|j = {(t, x) e Bfl(O) X (0, T); git, x) > e} . 

From Egorov's theorem and the compactness of g" and g" (Lemma 3.1), we have 
the existence of a set C with measure \ < on which and 
converge uniformly to g. Then, for n sufficiently large. 



g^ > e/2 in C, 
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and since 



□ 



(3.8) 



we can pass to the limit on C,, to deduce 

g ^ 

Since this holds for all 77 > 0, we can conclude 

m = uQy:, in 
for every R and e. We conclude that, 

m = ug^p on {g > 0} . 

Proof of Theorem 1.1: The weak formulation of (3.1) reads 
/ f"{A + V ■ Vx-iA - Vj;'J'V„'0) dvdxdt 

^E^d 

:= + 1^-1 ■) dvdx, e ((0, T) x R^^), 

where we have introduced the quantities 

/" = - / / ^{x-y)p{x,v)f^(y,w){w-v)S/ti'ip{x,v)dwdydvdxdt, 
Jo Jr*'^ 

I2 ^ - I I riu" - w)V„V dvdxdt. 

Jq JR2d 

By virtue of (3.2), we can pass to the limit in (3.8) to conclude 

rT 

fillet + V ■ '^xip - Vx'^Vyip) dvdxdt 

(3.9) 

h + hm r^' - / MiO, •) dvdx, 

where /i = - J^^^^ix ~ y)fix,v)f{y,w){w - v)Vvi^{x,v) dwdydvdxdt. 

From Lemma 3.2, we have that /"u" ^ fu in L°°((0, T); L«(M2'i)), for any 
q < , and hence there is no problems with passing to the limit in to discover 



lim I2 — — hm j f"-{u^ — v)V.t,ip dvdx — ~ / f{u — v)Vyil;dvdx. 



By setting this in (3.9) and recalling that ^I^ is arbitrary, we conclude that the limit 
/ is a weak solution to 

ft + div,(«/) - div,(/V,*) + div,(/F[/]) + div„(/(u - v)) = 0. 

This concludes the proof of Theorem 1.1. 



□ 
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